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In this paper, the gravitational search algorithm (GSA) is proposed as a 
method for controlling the opening and closing of airplane wing tires. The 
GSA is used to find the optimum proportional-integral-derivative (PID) 
controller, which controls the wing tires during take-off and landing. In 
addition, the GSA is suggested as an approach for overcoming the absence of 


the transfer function, which is usually required to design the optimum PID. 





The use of the GSA is expected to improve the system. Two of the most 
Keywords: popular optimisation algorithms-the harmony search algorithm (HSA) and the 
particle swarm optimisation (PSO)-were used for the sake of comparison. 
Moreover, the GSA-, HSA- and PSO-based optimum PID controllers were 
compared with one of the most important PID tuning methods, the Ziegler- 
Nichols (ZN) method. In this study, the integral time absolute error (ITAE) 
was used as a fitness function. First, four transfer functions for different 
applications were used to compare the performance of the GSA-based PID 
(PID-GSA), HSA-based PID (PID-HSA), PSO-based PID (PID-PSO) and 
Ziegler-Nichols-based PID (PID-ZN). Next, the GSA was used to design the 
optimum PID controller for the opening and closing systems of the airplane 
wing tires. The results reveal that the GSA provides better outcomes in terms 
of ITAE when compared with the other adopted algorithms. Furthermore, the 
GSA demonstrates a fast and robust response to reference variation. 
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1. INTRODUCTION 

Controlling the wing tires is an important aspect of airplane take-off and landing. There are many 
challenges associated with the control process, such as settling time, rise time, overshoot and steady-state 
error [1]. The goal is always to enhance the transient and steady-state performance of a system [2]. 
Accordingly, the design of a robust controller for the purpose of opening and closing the wing tires is an 
important matter [3]. One of the most important controllers reported in the literature is the PID controller [4]. 
The performance of the PID controller depends on three parameters: proportional gain (K_p), integral gain 
(K_i) and derivative gain (K_d) [5]. These parameters produce an expected transient response which 
corresponds to each operating point [6]. A first and integral step in designing the PID controller is deriving 
an appropriate mathematical model; after this, various techniques can be implemented to determine PID 
parameters (K_p, K_i and K_d). In the literature, the ZN technique is the one of the most important for 
designing the PID controller [7]. 
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However, the derivation process for the mathematical model is complicated-sometimes impossible. 
Recently, however, complex, real-world problems have been solved using computational intelligence 
optimisation algorithms (CIOAs) [8]. The term ‘computational intelligence’ indicates solving a problem 
using artificial or natural populations, and suggests that complex problems can be solved using CIOAs. Thus, 
CIOAs are suggested in this paper in the absence of a mathematical model. Furthermore, the use of CIOAs is 
expected to improve the opening and closing system of the wing tires. A CIOA can often solve a problem 
more quickly, or find an approximate solution when classical methods fail, especially in the case of 
multimodal optimisation problems [9]. Many CIOAs are mentioned in the literature. These algorithms can be 
divided in two groups: evolutionary and swarm intelligence algorithms [10]. 

Darwin’s theory of evolution is fundamental to the evolutionary algorithmic group, the most popular 
paradigms of which are the genetic algorithm (GA) [11], evolutionary programming [12], differential 
evolution [13], evolutionary strategy [14] and genetic programming [15]. However, falling into local minima 
and difficulty finding the optimum solution to the multimodal problem are the main drawbacks of these 
algorithms. 

The second group is swarm intelligence, the algorithms of which mimic the social behaviour of 
animals and insects [16]. PSO [17], ant colony optimisation (ACO) [18] and artificial bee colony optimisation 
(ABC) [19] are the main paradigms. The PSO algorithm is one of the modern heuristic algorithms, it was 
enlivened by the social behavior of bird and fish tutoring and has been discovered to be vigorous in solving 
continuous nonlinear optimization problems [20]. Where, PSO algorithm is applied to improve the surface 
grinding process parameter in both rough and finish grinding condition. The The ideal worth of the decision 
factors were acquired by utilizing the mathematical model with the purpose of decreasing the production cost 
and expanding the production rate as well as improving the surface finish [21]. 

The main limiting factors of these algorithms are early convergence and stagnation in local minima. 
Other swam intelligence optimisation algorithms have been developed such as such as HSA [22], GSA [23], 
firefly algorithm (FFA) [24] and bat algorithm (BA) [25]. This study aims to present a method for finding the 
optimum PID controller in the cases of both absence and presence the transfer function of any system. For 
this purpose, a well-known meta-heuristic optimisation called GSA is used. This method aims to improve the 
performance of any system. Furthermore, it is intended to solve the problem of trial and error in determining 
the parameters of the PID controller, which are K_p, K_i and K_d. Therefore, the use of the GSA is expected 
to enhance the system response. For this purpose, a formulated fitness function, ITAE, is used. The GSA 
aims to reduce ITAE by regulating the control parameters (K_p, K_i andK_d) of PID controller during the 
design phase. 


2. PROPOSED CONTROLLER OVER WINGS TIRES 

The most important element of the control process of the wing tires is the proper selection of the 
PID parameters. The improper selection of these parameters will lead to a bad response in the process of 
opening and closing the wing tires. Accordingly, this paper presents a methodology for designing the PID 
controller without a mathematical model, instead using meta-heuristic optimisation algorithms. This process 
will avoid the trial and error process in the design of the controller. Moreover, this methodology is expected 
to enhance the performance of the opening and closing process of the wing tires. 


2.1. Proposed meta-heuristic optimisation algorithms 

In recent years, meta-heuristic optimisation algorithms have been used to solve problems in various 
fields, such as commerce, economics and engineering. Meta-heuristic optimisation algorithms are often used 
based on their ability to address complex problems in real-world situations; ease of use and efficiency also 
contribute to their attractiveness. In this work, a GSA has been developed to enhance the opening and closing 
process in wing tires. Through the use of this algorithm, an improvement in the performance of the airplane 
wing tires is expected. 


2.2. Formulation of the optimal PID controller 

The formulation process for the optimisation-based wing tires consisted of three parts: the 
determination of the search space size; the formulation of a suitable fitness function for the system; and the 
identification of optimisation limitations. The optimisation algorithm search the optimum value of the search 
space through depending on the fitness function and optimisation controls. 


2.2.1. Size of the search space 


The size of the search space means the input variables through which the optimisation algorithm 
should find their optimal solutions. The size of the search space can be expressed as: 
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lj = [KË ik | (1) 


The J; j is the jth solution in the agents through the ith iteration, ki J is the tth element of Ii ;, and n is the total 


number of variables. 


j? 


2.2.2. Fitness function 

The second component of the optimisation process is the fitness function, which is important for 
evaluating the performance of the optimisation algorithm. Accordingly, the fitness function has been 
developed to find the optimum design of the PID controller in such a way that the I; ; generates the optimum 
control action. ITAE has been developed as the fitness function, and is found by taking the difference 
between the reference input (Yin) and the measured output (Yput). The ITAE is expressed by: 


ITAE = f° t\Vin(t) — You (t) lat (2) 


Where Yin is the reference input, Y,,,- is the measured output and s is the sample number. 
During the optimisation process, a minimal fitness function indicates an effective system; thus, the 
optimisation algorithm is aiming to minimize ITAE. 


3. PROPOSED CONTROLLER-BASED WING TIRES 

To demonstrate the effect of the proposed controller on the opening and closing of the airplane wing tires, 
the system has been modelled using MATLAB environment, as shown in Figure 1. The output signal (Yput) is 
measured and compared with the input signal (Yin). The difference signal between Yin and Yjy¢ is fed into the 
optimisation algorithm to find the optimumK,, K; and K4 of the PID controller. As mentioned before, the GSA has 
been used to enhance the opening and closing of the wing tires. The performance of the GSA has been compared 
with HAS and PSO. Each algorithm has special variables that differ from those of other algorithms. These 
parameters must be settled in the initialisation stage. The algorithm parameters are listed in Table 1. 


Optimization 
Algorithm 





Wing Tires 





Figure 1. Closed-loop system of the wing tires 


Table 1. Parameters of the optimisation algorithms 








Parameter GSA HAS PSO 
Agent Size (N) 20 20 20 
Maximum Iteration (T) 500 500 500 
Initial gravity 100 - - 
E 2.2204e-16 - - 
Harmony consideration rate - 0.9 
Minimum pitch adjusting rate - 0.4 
Maximum pitch adjusting rate - 0.9 
Minimum bandwidth - 0.0001 
Maximum bandwidth - 1 - 
ai , Pf - - 2.5 
Bi, ap - - 0.5 





4. IMPLEMENTATION OF THE GSA 
As mentioned before, the GSA has been used to improve the opening and closing system of an 
airplane’s wing tires. Introduced by Rashedi et al. [23], the GSA is a meta-heuristic optimisation algorithm 
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inspired by the law of gravity. Accordingly, this algorithm is based on Newtonian gravity. The 
implementation of the GSA in the design process for the PID controller-based wing tires begins by resetting 
the GSA parameters. These parameters are number of agent (N), number of iteration (T), initial gravity and 
final Kbest. The initial agents are encoded using (1). The evaluations of these agents are achieved using (3). 
Each agent is updated through the (3): 

Xow = xf + View (3) 
Where velocity View) of each agent is calculated as: 

View = rand; x Vf + af (4) 


The acceleration of the ith agent and dth direction (af) is computed as: 
Fi 
af = (5) 


Where M; represents the inertia of the ith agent. The stochastic characteristic of the GSA algorithm has been 
achieved using the total force. The total force of the ith agent and dth direction (FÊ) is represented as: 


Ff = Èj- jei rand Fg (6) 


The force acting of the ith mass from jth mass (Fi) is defined as: 





d — 7 MpixMaj( yd _ yd 
MaG errs S (7) 
Where G, Mpi Maj, and € are the gravitational constant, passive gravitational mass of the ith agent and active 


gravitational mass of the jth agent, respectively. The euclidian distance between ith agent and jth agent (R;;) 
is defined as: 


Ry = XE, XL, 8) 


This process is repeated until the maximum iteration is reached, as shown in Figure 2 in appendix. 


5. RESULTS AND DISCUSSION 

Two tests have been carried out in this study. Presented in section 5.1, the first test shows the 
performance of the GSA in designing the PID controller for some plants. The second test, presented in 
section 5.2, shows the performance of the GSA in solving the problem of trial and error procedure when 
designing the PID controller for opening and closing systems for airplane wing tires. 


5.1. Test-1 

The GSA has been used to design an optimum PID controller. The GSA has also been compared 
with HSA and PSO. For further comparison, the ZN method, a renowned PID tuning method, has also been 
used in the design of optimum PID controllers. Four different transfer functions, shown in Table 2, have been 
used to evaluate the four controllers (PID-GSA, PID-HSA, PID-PSO, and PID-ZN). The results of these 
controllers, based on ITAE, are shown in Table 3, and the best performance is bolded. Table 3 clearly 
indicates that the GSA, HAS and PSO have achieved minimum ITAE when compared to ZN. Furthermore, 
Table 3 shows that the performance of the PID-GSA is better than the other controller for all transfer 
functions. This highlights the robustness of the PID-GSA. 

The convergence characteristics of the GSA-, HAS- and PSO-based first-transfer function (G,) are 
shown in Figure 3. Compared to the other algorithms, the GSA achieved fast convergence. Thus, the GSA 
has superior convergence characteristics for this transfer function. Furthermore, the GSA achieves minimum 
ITAE compared to the HSA and PSO. The step response of the G4 is shown in Figure 4, and shows that the 
PID-GSA first-transfer function performs better than the other controllers. The PID-GSA achieves faster rise 
time, minimum overshoot and minimum peak time, indicating high efficiency. 

The comparison between PID-GSA, PID-HSA, PID-PSO and PID-ZN for four transfer functions is 
presented in Table 4. In a control system, rise time, settling time, overshoot and peak time are the main 
parameters considered; these are shown in Table 4. The best performances are bolded. For all tested transfer 
functions, the GSA achieved the best results. This underscores the validity of the GSA and its ability to deal 
with different types of transfer function. 
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Table 2. Tested transfer functions 


ID Plant Name Transfer function of the plant 
0.2304s + 0.01265 














1 G(s) Heat exchanger Temperature G,(s) = =a 0.066895 T 0005744 
P i G. = : 
2 G(s) Satellite model 2(s) (5 + 2.723)(s + 2432) 
3 G3(s) Water tub G3( ja 0 
as E 3S. = 3s + 1)(29s +1) 
5.54 
4 G,(s) Potential hydrogen neutralization G,(s) = 220 F1 





Table 3. Results of all optimisation algorithms for all plants 

ID Plant Method Kp Ki Kg ITAE 

1 Heat exchanger PID-GSA 22.3663 9.0319 1.0912 1.6990 
PID-HSA 13.6852 9.9338 1.0180 1.8923 
PID-PSO 12.0639 9.8883 1.0203 1.9740 
PID-ZN 2.8272 1 0.25 13.8320 

2 Satellite model PID-GSA 32.8582 22.9864 3.6112 0.4774 
PID-HSA 21.8426 19.6855 2.3160 0.6243 
PID-PSO 17.8358 18.6048 1.5826 0.7028 
PID-ZN 30 20 5 0.9989 

3 Water tube PID-GSA 157.5407 41.2388 3.9137 5.8112 
PID-HSA 105.2762 1.6099 28.9902 6.2834 
PID-PSO 100.3626 31.32554 11.25852 7.0740 








PID-ZN 78 1.5 0.375 23.2565 
4 Potential PID-GSA 19.6946 1.21654 2.65325 0.5864 
hydrogen PID-HSA 9.3999 0.8712 1.6747 0.8425 


neutralization PID-PSO 19.4398 1.1738 2.1015 0.9495 
PID-ZN 1.05880 1.3240 0.2120 3.9760 











1.4F F T T 

















2 
£ 
[e] 
S 
£ 
T 
a 
15! : : : i : ; 
0 100 200 300 400 500 10 
Iteration Time (s) 
Figure 3. Convergence characteristics of the heat Figure 4. Step response for the heat exchanger 
exchanger transfer function (G1) transfer function (G1) 


Table 4. Comparison between responses of the plants based on different methods 








ID Plant Method Rise time Settling time Over shoot Peak time 
1 Heat exchanger PID-GSA 0.49 s 4.2 s 6.6% 1.06 s 
PID-HSA 0.64 s 3.9 s 14.1% 1.53 s 
PID-PSO 0.69 s 3.8 s 16.5% 1.69 s 
PID-ZN 2.07 s 20s 20% 4.81s 
2 Satellite model PID-GSA 0.41 s 0.9 s 2.1% 0.68 s 
PID-HSA 0.49 s 1.3 s 3.3% 0.84 s 
PID-PSO 0.53 s 14s 4.9% 0.92 s 
PID-ZN 0.56 s 3.55 1.4% 15s 
3 Water tube PID-GSA 1.24 s 13.24 s 36.8% 2.26 s 
PID-HSA 1.53s 18.72 s 35.5% 2.74 s 
PID-PSO 1.53 s 22.73 s 43.9% 2.87 s 
PID-ZN 1.61 s 26.63 s 56.8% 3.12 s 
4 Potential Hydrogen PID-GSA 0.3 s 14s 16% 0.5 s 
neutralization PID-HSA 0.42 s 1.9 s 18.7% 0.72 s 
PID-PSO 0.31 s 17s 24.5% 0.5 s 
PID-ZN 0.75 3.758 16.1% 1.6s 
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5.2. Test-2 

In the second test, the GSA is used to control the process of opening and closing the airplane wing 
tires. The convergence characteristics of the GSA-, HAS- and PSO-based wing tire problems are exhibited in 
Figure 5. The figure shows that the GSA provides better solution than the other algorithms. 

The performance of the GSA was compared with HAS and PSO, as well as optimum PID through 
MATLAB. Four controllers, PID-GSA, PID-HSA, PID-PSO and PID, are illustrated to deal with opening 
and closing of the wing tires of the airplane. The first test is to change the wing tire position from close to 
open to prepare for the landing as shown in Figure 6. The performance of the four (PID-GSA, PID-has, PID- 
PSO and PID) controller-based landing situations are shown in Figure 7, which clearly indicates the best 
performance was achieved by PID-GSA. The PID-GSA succeeds in achieving faster steady-state error, rise 
time and settling time than the other controllers. 
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Figure 5. Convergence characteristics of Figure 6. Change in wing tire 
GSA-HAS-and PSO-based wing tire problems position from, (a) close to (b) open 
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Figure 7. Response of wing tires in landing situation (from 90° to 0°) 


The performance comparison of the PID-GSA, PID-HSA, PID-PSO and PID controller-based wing 
tire problems is shown in Table 5. Again, ITAE is used as the fitness function to evaluate all the controllers. 
The best performance is bolded. Table 5 shows that the GSA succeeds in providing better PID controller 
where minimum ITAE has been achieved. 


Table 5. Results of four controller-based wing tire problems 
Method Ky K; Kg ITAE 
PID-GSA 2.999345097233e+05 1.159227855138  4.097434603125e+04 1.1830e+03 
PID-HSA 2.610763651923e+05 1.863932871596 — 5.653000411136e+04 1.194e+03 
PID-PSO 2.754962374074e+05 2.096996978282  4.764440518909e+04 1.230e+03 
PID 2e5 10 5e4 3.2375e+03 
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6. CONCLUSION 

Meta-heuristic optimisation algorithms have been proposed in many applications to solve complex 
problems faster and more effectively. This study introduces a method for using a meta-heuristic optimisation 
algorithm, GSA, in the design of an optimum PID to open and close airplane wing tires. The meta-heuristic 
optimisation algorithm is applied to solve the problem in the absence of a mathematical model. GSA use also 
leads to improved system performance. In the study, the GSA is compared to both HAS and PSO, the PID- 
GSA, PID-HAS and PID-PSO, as well as the PID-ZN. Two tested were executed. The first test uses four 
transfer functions, while the second uses the GSA, HAS and PSO to design the optimum PID controller for 
the opening and closing system of the airplane wing tires. The ITAE has been formulated as a fitness 
function; the meta-heuristic optimisation algorithms aim to reduce ITAE during the design phase. The 
PID-GSA achieves better performance based on all transfer functions. Moreover, the PID-GSA achieves 
robust and favourable results for opening and closing the airplane wing tires of. Therefore, the PID-GSA can 
be effectively applied for use with the airplane wing tires. 


APPENDIX 


Generate initial population using (1), Agent Size 


(N), Maximum Iteration (T) 

















No 


Run the system 


Calculate objective function (ITAE) using (2) 





Calculate the distance using (9) 
Calculate the total force and force using (7,8) 


Calculate the acceleration using (6) 


Calculate the velocity using (5) 
Update the agents using (4) 


Reach Maximum population? 


Calculate objective function (ITAE) using (2) 


Objective function evaluation 


Output the optimal parameters 


Figure 2. Implementation of the proposed GSA 
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